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ABSTRACT 

For cutoff potentials, a condition which is not a limitation for the calculation of physical 
systems, the S-matrix is meromorphic. We can express it in terms of its poles, and 
then calculate the quantum mechanical second virial coefficient of a neutral gas. 

Here, we take another look at this approach, and discuss the feasibility, attraction 
and problems of the method. Among concerns are the rate of convergence of the 'pole' 
expansion and the physical significance of the 'higher' poles. 



1 Introduction 



This work has had a long gestation. Our earliest notes date from July 1976! At that 
time, two of us (M. Berrondo and S.Y. Larsen) obtained what we thought was an 
attractive formula for the q.m. second virial coefficient in terms of the location of 
its poles, as an alternative to its well known M and customary formulation in terms 
of phase shifts and bound state energies. Soon after doing this work, and obtaining 
some results, we were made aware that H. Nussenzveig had already published^ an 
attractive article on the same subject and we stopped our efforts. 

Recently, however, there has been interest in some other quarters0, related to the 
pole expansion of the S-matrix, and we thought that useful details of such expansions 
might be properly be brought to light within the context of the statistical mechanics of 
the virial. Together with other colleagues, we look at the rate of convergence for phase 
shifts and for the virial. We look at asymptotic expansions for the poles. We try to 
see whether we can use tricks to either accelerate the convergence or sum background 
poles. We look at different model potentials. In the case of hard spheres, we obtain an 
expansion for the virial valid at low temperature. 



2 The second virial in terms of the poles 

For potentials decreasing faster than exponentials, like gaussian or cutoff potentials 
the Jost function is an entire function and, accordingly, the S matrix is meromorphic. 
We can express it as a exponential together with a product of poles : 

J2iS, i \£ -2ika tt kn,£ + k /1A 



S e (k) = e 2ibl = (-Y e~ 2ika J] 



k n i k 



where a denotes the range of the potential and the k n /s are the poles labelled in order 
of increasing modulus. It should be noted that this expression, apart from the factor 
which we will examine later, is equivalent to the one derived some years ago by 
Nussenzveig @] 

St(k) = e ~ 2 ^ II ) P Z \ ■ ( 2 ) 

We recall that the original quantum mechanical formulation of the Boltzmann part 
of second virial coefficient, in term of phase shifts and bound state energies, reads 



(£2 W = -2^\%N 



V(2£+l)e~ 2/3 ^ + - r dke~ 2l5k2 V(2£ + 1)4W 
bp 7r Jo , dk 



(3) 

where f3 = 1/(kT) in terms of the Boltzmann 's constant k and the temperature T, 
At denotes h/ V 2itmKT and J\f is the number of particles in the volume V. (Note that 
2/3 = A T /(27r)). The sum for the bound states runs over £ and the number of bound 
states for each £. 
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In later calculations and formulations |J, a partial integration has often been per- 
formed to yield the virial in terms of the phase shifts, themselves, and the bound state 
energies. Also for convenience, we have divided the virial into two parts, a Boltzmann 
plus an exchange part. We focus on the Boltzmann part, but can at any moment obtain 
the exchange part by minor internal changes of sign and the inclusion of a perfect gas 
term. 

Since the S matrix, here, has a compact expression in terms of its poles the previous 
equation also has such an expression. Indeed the derivative d5e/dk is nothing else than 



1 



k n ,t ~\~ k k n ^i k 
Introducing (4) in (3) we look at the contribution from the bracket 

k n ,£ 



(4) 



£(2£ + IK 2 ^ + 1 £(2£ + 1) (~J^ - E ~ f°° dk ^ 



U2 _ U2 



(5) 



The integral involving the pole expansion can be written in terms of the error function 

2 



erf(z) = -= f e~ v2 dv . (6) 

v7T JO 



In Appendix A we show that 

- ¥ f * e ~ m * whi, = r ( erf (-*».«v^) =f i) . m 

with the sign - when the $s(k n ,e) is positive and the sign + when the $s(k n ,e) is negative. 
The bracket then reads 

J2(2i + l)e~ 2 ^ + £(2£ + 1) (~^= + \ E U{-ik n ^) T 1 

B,£ i \ v2V4vr^ l n L 

(8) 

We now put together the terms corresponding to bound states (52b i) an d the terms 
corresponding to bound states, in the expression involving the poles (J2 n ,e) • The poles 
involved in these terms are situated in the upper half plane $s(k) > 0. According to 
our previous discussion, the exponential exp(— 2[3{k n j) 2 ) is then weighted by the factor 
(-1/2). Combined with the exponential terms of the bound states, this yields the factor 
+ 1/2 as for the other poles. The bracket then reads 

Remembering, now, that erf(— z) = — erf(z) and 1 — erf (z) = erfc(z), and using 2(3 = 
\t/(2tt) we obtain our final expression ( in terms of At): 



{B,) Boltz = -2^ T N E(2^ + l^ 



At 1 , o s „ , . At 



-^exp(--5/c 2 i£ )erfc( 



V2X T ' 2^^ K 27t' n ' a K V2^ 



(10) 
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3 The phase shifts in terms of poles 



Clearly our procedure becomes more attractive if few poles are required to reproduce 
the second virial to good accuracy. In this section, we examine how the expansion, 
Eq.(4), reproduces the phase-shifts. 

For a pure hard sphere potential, the answer is very pleasing: the number of poles 
for the £ th partial wave equals £. 

For £ = 0, the phase shift of a hard sphere of radius a is just — ka and the poles 
do not contribute. For the higher ts, we find that our formula for the S-matrix works 
perfectly, but that the factor (— l) e is not needed. We note that asymptotically the 
phases tend to — ka — in/2. 

Unfortunately, for the more common cutoff potentials, the number of poles appear- 
ing in the expansion is infinite and, practically, to obtain the phase shift to, say 5 
digits accuracy, the number of poles required is impressively large. We note, though, 
that except for poles found on the imaginary axis, they occur in pairs, in the third 
and fourth quadrant, and thus for these pairs, it is sufficient to determine the poles in 
the last quadrant. For a potential made up of a (repulsive) hard core plus an attached 
attractive square well, we find the results presented in table 1. 



N 


5o 


So 


22 


.067181 


(.9909869, .1339585) 


72 


.073081 


(.9893373, .1456426) 


122 


.074379 


(.9889561, .1482090) 


522 


.076023 


(.9884633, .1514605) 


1022 


.076312 


(.9883756, .1520316) 


2022 


.076473 


(.9883265, .1523510) 


10022 


.076620 


(.9882818, .1526403) 


50022 


.076654 


(.9882713, .1527083) 


150022 


.076661 


(.9882693, .1527210) 


Exact 


.076664 


(.9882682, .1527283) 



Table 1. Poles for 1 antibound state + N pairs 



This is for t = 0, ka = q = 0.1, the de Boer parameter A* = [h 2 /mVoa 2 ) 1 ^ 2 equal 
to 10, the hard core radius a and finally a, the outer limit of the attractive potential, 
equal to 2.85o\ 

The expansion converges but so slowly that we are faced with the necessity of 
increasing the rate of convergence of the series. We note that at the origin (k — 0), 
the phase shift behaves like k 2l+1 . Thus, for a given £, derivatives up to 2£ are equal 
to zero. We obtain: 

- a +-Y,T~ = for £>1 

E72J+1 =° l<j<t-lfoi t>2. (11) 

n k n e 
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We can now try to accelerate the convergence of (4) by subtracting the terms above, 
multiplied by appropriate powers of k. We obtain, for t > 1, 



dk 1 



which can also be written as 

d_ 

dk' 



+ 



£-1 

E 



k 



2.7 ' 



k n ,e + k k nt f k k n ^ - y k n 



(12) 




U2 _ h.2 

K n,t ft 



(13) 



We see that, for high orders, the terms behave as 1 / (k nt p) 2£+z . 
Integrating the above equations with respect to k we find 



In 



A: 



k n ,e + k i x - - / n 
^k) ~ f^ Q 27+1 ^ 



A: 



> 1 



(14) 



For £ = l,q — 0.1 and 23 pairs of poles, we obtain 1.478248 10 3 compared to the 
exact result of 1.47826613 10~ 3 and, similarly, for I = 2, 4.04298921 10" 6 instead of 
4.042988889 10~ 6 . The method, however, deteriorates as the energy increases. 

To remedy the slow convergence for i = 0, we introduce the derivative of the phase 
shift (non zero) at k — , i.e. 



+ 



k n ,o + k k n fl k k n fl 



(15) 



We then easily calculate the new term using trigonometric functions. 

To prevent the deterioration as k increases, we limit the number of subtractions to 
2. The results are then presented in table 2. Clearly, the number of poles required, 
I being fixed, rises as the energy increases. On the other hand, at fixed energy, the 
number of poles needed is least when I is largest. 



ka I 





1 


2 


3 


4 


5 


6 


1 


69 


70 


19 


20 


21 


22 


23 


2 


141 


142 


37 


38 


39 


40 


41 


3 


217 


218 


55 


56 


57 


58 


59 


4 


291 


292 


73 


74 


75 


76 


77 


5 


367 


368 


91 


92 


93 


94 


95 


6 


443 


444 


111 


112 


113 


114 


115 


7 


519 


520 


129 


130 


131 


132 


133 


8 


597 


598 


147 


148 


149 


150 


151 


9 


673 


674 


165 


166 


167 


168 


169 


10 


751 


752 


185 


186 


187 


188 


189 



Table 2. Poles required for 5 digit accuracy 
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4 Asymptotic formulas 



The simplest way to remedy a lack (or a limited number) of poles is to determine an 
asymptotic formula. Nussenzveig, in his bookQ], derived such an asymptotic expres- 
sion, which we extended to 'handle ' hard cores, but we need to increase its accuracy. 

We begin by looking for an asymptotic formula for the s-wave, and then generalize 
to higher angular momenta. The zeros of the Jost function F(k) are the poles of the 
S-matrix F(—k)/F(k). They are given (see Appendix B) by solving the equation 



D 2ik(a—o) 



+ 



Ak 2 

W) 

1 



1 + — \2M 
2k \ 



V'{a) 
V(a) 



'V"(a) 



Ak 2 \ V(a) V(a) 



Via) 2 , . V'(a) 9 N 

- 2V(a) + 2M^4- - 2M 2 



V(a) 



(16) 



where M denotes M — J° V(r') dr' (for a potential finite at the origin put a = 0) and 
V(a) is assumed different from zero. 

We thus obtain some corrections to Nussenzveig's leading term, which were calculated 
from the equation: 

Ak 2 



3 2ifc(a— a) 



V(a) 



(17) 



To get the correct value we first solve equation (16) by setting 

k (a — a) = rnr — e n/2 — i A , (11 
where e = 0, 1 according to whether V(a) is positive or negative, and, then iterating 

a _ 2 0r o - iA) 



A 



with x = Tin — e n/2 and A 2 = \V(a)\(a — a) 2 . The first correction is 

A = In — , 



(19) 



(20) 



which corresponds to that of Nussenzveig. The third iteration, practically, provides us 
with the exact value so that 



A ~ A 2 = In [(2x - 2iln((2x - 2iln(2x /A))/A)) /A] . 



(21) 



To get higher order corrections we have to include additional terms in (17), involving 
higher powers of 1/k, and solve by iteration 



A _ xp - iA 
6 A 



1 + 



ax 



+ 



Oi 2 



with 



Oil 



a 2 



i (a — at 



[a — a) 



2M 



x — iA (xq — iA) 2 
V(a) 



+ 



(22) 



V"{a) 3 V'(af 



V(a) 4 V(af 



2V(a) + M 



V(a) 



M 



(23) 
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To test the method for the potential that we used before ( hard core + square well), 
we define A = 2-n-(a — a)/A*. An expansion in 1/ A; up to the 4 th order yields the results 
in table 3. We see that we are able to reproduce the exact results, for the higher poles, 
to 5 digits. 





n 


k n ,0 


appr. 

exact 


20 
20 


(33.0670568,-2.5190331) 
(33.0670586,-2.5190330) 


appr. 

exact 


21 
21 


(34.7670884,-2.5459707) 
(34.7670900,-2.5459707) 


appr. 

exact 


22 
22 


(36.4669649,-2.5716312) 
(36.4669647,-2.5716312) 


appr. 

exact 


23 
23 


(38.1667055,-2.5961305) 
(38.1667061,-2.5961304) 



Table 3. A* = 10, a/a = 2.85 



For waves with higher values of £, we show, in Appendix C, that we can generalize the 
approach of Appendix B, for example by setting the free Jost solution equal to 

We(kr)=^ 1 ] f^rH^ 1/2 (kr), (24) 

and proceeding in a manner similar to that used for £ = 0. 

In table 4, we examine the quality of our approximations for values of i not equal to 
zero. Our results for t = 1,2 are of comparable quality to those that we found for the 
I = 0. The key parameter is the ratio tjkia — a). So long as it is small the asymptotic 
results will be good. We see that for large ts, we have to proceed to a larger value of 
n, before the behaviour of the poles becomes asymptotic. 



i 




^23/ 


1 


appr. 

exact 


(38.1758622,-2.5955925) 
(38.1758614,-2.5955925) 


2 


appr. 

exact 


(36.4956875,-2.5698867) 
(36.4957008,-2.5698843) 


3 


appr. 

exact 


(36.5243519,-2.56815261) 
(36.5244255,-2.56813956) 


4 


appr. 

exact 


(34.8672057,-2.53969925) 
(34.8674660,-2.53965187) 


5 


appr. 

exact 


(34.9169320,-2.5366040) 
(34.9175873,-2.5365038) 


6 


appr. 

exact 


(33.2864742,-2.5046694) 
(33.2883530,-2.5045969) 



Table 4. A* = 10, a/a = 2.85 
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When the potential has no hard core, we show, in Appendix C, that the asymptotic 
expression of the poles is obtained (in lowest order) from the solution of 



e , ika = { _ Y m . (25) 

We then recover the (— ) e dependence, mentioned earlier by NussenzveigQ]. This de- 
pendence disappears when the potential incorporates a hard core. 



5 Virial 

Noting, in the previous section, the necessity of accelerating the convergence of the 
pole expansion, we present here the formalism for doing this for the virial. Afterwards 
we will discuss the application to various potentials. 

Given the slowness of the basic expansion in terms of poles, we modify the basic 
virial equations, proposing two different versions. In the first one we simply write 



k=0 



1 yj — 

i n k n ,£ 



and add and subtract this from the derivative expression. This yields: 

1 



(B 2 ) 



Boltz 



+ 2^+1) (' 



2 1 / 2 X T 
exp(- 



k=0 



X T U 2 ^ _f„/„- X T , 

-At, 



2tt 



K,i) erfc O 



2,71 



J + 



2 1 / 2 i ' 



We have used the fact that 



' d_ 
dk 



S t (k) 



for e > . 



k=0 



We can push this further, using Eq.(ll) of our paper. We then find that 

1 



B 2 



Boltz 



-2 1 ' 2 \%M 



2V2A 



T 



k=0 



+ \ £(2£ + 1) exp(-^ ) erfc(i jW 
2 . 2tt V2tt 



+ 



1 



e-i 



2V2A 



E( 2 ^ + !)— + E( 2 ^ + 1) E E " !) 



7T 



(26) 



(27) 



(2* 



(29) 



' \*.„ '-I- n j=l A T fcj^ 1 



Finally, what we calculate is this virial divided by that obtained classically for a pure 
hard core of radius a, i.e. 



[Bl) 



Boltz 



27TA/V 3 



(B 2 ) 



Boltz 



(30) 
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5.1 Hard spheres 

We illustrate in table 5, how with a limited number of poles (about 20) we reproduce 
the 2 nd virial coefficient for hard spheres, over a range of relatively low temperatures, 
for which these formalisms are designed. 





Xt/ct = 1 






Xt/ct 


= 30 


t 


BLK 


1 




Poles 


BLK 


I 




Poles 





.4774648 


29276 


.477464829276 


.42971834634 


8E+03 


.429718346348E+03 


1 


.75234870^ 


^365 


.752348708365 


.445367908^ 


373E+01 


.445367908873E+01 


2 


.584467004923 


.584467004923 


.8675880964^ 


^8E-02 


.867588096488E-02 


3 


.3359358 


32978 


.335935832978 


.848769265246E-05 


.848769265246E-05 


4 


.1554667 


87903 


.155466787903 


.544399040752E-08 


.544399040752E-08 


5 


.059852912973 


.059852912973 


.258124235604E-11 


.258124235604E-11 



Table 5. The hard sphere virial, from Eq.(29) 



In Appendix D, we outline how we recover terms in a low temperature expansion 
(ct/At small), a number of which are found in BLK. 

5.2 Hard sphere + square well 

For A* = 4, the 2 particles have a bound state for I = and one for i — 1. These 
bound states then dominate the very low temperature behaviour of these partial virials 
(i.e. for these angular momenta), and therefore for these virials. This is due to the 
large exponential terms that arise for large Xt/ct. This already becomes evident for 
Xt/ct = 2, and becomes more important for Xt/ct = 5, 10, 20, 30 • • • 

Still, for Xt/ct = 2, it is useful to see the virial calculated and displayed as function 
of the number of poles taken into account, for different values of I. In table 6 we 
subtract the first derivative as in Eq. (27). 

Here it is well to comment on our reference BLK. The authors BLK published 
results on hard spheres||, many years ago, but not partial results for distinct values 
of £, nor results for a hard sphere + a square well. Here, by BLK, we mean that one 
of the previous authors (SYL) is using some of the old programs and the old methods 
to obtain results which can be used to benchmark the use of the poles. These results, 
however, have not been obtained with the same need and desire for accuracy that the 
old work required. The programs have now been used with machines with smaller word 
lengths, the calculations done with larger meshes, etc. The values of the bound states 
have also been obtained from our work with the poles. 
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N i 





1 


2 


3 


B+AB 


-4.96245020 


-13.49991189 


9.568499294 




1 


-5.38351004 


-14.59384163 


-13.71610648 


7.071340875 


2 


-5.42559329 


-14.70879231 


-15.01832365 


-8.916442069 


5 


-5.44089652 


-14.75294247 


-15.23836884 


-10.22935288 


10 


-5.44304631 


-14.75932474 


-15.25610946 


-10.27349937 


20 


-5.44347274 


-14.76060000 


-15.25890823 


-10.27875744 


50 


-5.44355258 


-14.76083933 


-15.25936689 


-10.27949972 


100 


-5.44355890 


-14.76085827 


-15.25940041 


-10.27954958 


200 


-5.44355982 


-14.76086106 


-15.25940520 


-10.27955649 


300 


-5.44355993 


-14.76086137 


-15.25940573 


-10.27955725 


400 


-5.44355996 


-14.76086145 


-15.25940587 


-10.27955744 


BLK 


-5.44355995 


-14.7608616 


-15.2594062 


-10.2795575 




Table 6. 


A* = 4, a/a = 


= 2.85, \ T /a = 


2. 



The virial as a function of the number N of pairs, for different angular momenta. 
The contribution of bound (B) and antibound (AB) states is given separately at the 

beginning of the table. 

Here it is well to comment on our reference BLK. The authors BLK published 
results on hard spheres |]], many years ago, but not partial results for distinct values 
of £, nor results for a hard sphere + a square well. Here, by BLK, we mean that one 
of the previous authors (SYL) is using some of the old programs and the old methods 
to obtain results which can be used to benchmark the use of the poles. These results, 
however, have not been obtained with the same need and desire for accuracy that the 
old work required. The programs have now been used with machines with smaller word 
lengths, the calculations done with larger meshes, etc. The values of the bound states 
have also been obtained from our work with the poles. 

This said, we present two tables, exclusively with poles, for A* = 10 and a/a = 2.85, 
for which there are no bound states. In table 7 we show, for i = 0, the convergence of 
the virial as a function of the number of pairs, for different values of Xt/ct. 



N \ T /a 


1 


2 


5 


10 


AB 


0.598187603 


1.554449992 


0.649915614 


-20.719592323 


1 


0.503078797 


1.357448734 


0.454501918 


-20.909727874 


2 


0.457717502 


1.318856199 


0.420668117 


-20.942948517 


5 


0.434770286 


1.301162611 


0.404123797 


-20.959343787 


10 


0.431866213 


1.298475743 


0.401489560 


-20.961970746 


20 


0.431323513 


1.297945047 


0.400962083 


-20.962497768 


50 


0.431225758 


1.297847856 


0.400865048 


-20.962594781 


100 


0.431218271 


1.297840377 


0.400857571 


-20.962602258 


200 


0.431217192 


1.297839297 


0.400856492 


-20.962603337 


300 


0.431217071 


1.297839177 


0.400856371 


-20.962603457 


400 


0.431217040 


1.297839146 


0.400856349 


-20.962603488 



Table 7. A* = 10, a/a = 2.85, t = 
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In the table 8, we sum over the angular momenta, up to I = 7, noting that the 
largest contribution for I = 7 is of the order of 10~ 7 , and show the convergence of the 
virial as a function of the number of pairs per £, for different temperatures. 



N At/ct 


1 


2 


5 


10 


AB 


4.19599724 


9.11116433 


13 


85909885 


-5.12661171 


1 


10 


80533807 


15.47507182 


8 


97472495 


-18 


.26521543 


2 


11 


51743916 


10.90069896 


-0 


65662473 


-28 


.69768995 


5 


2 


21953160 


0.10665942 


-10 


39425326 


-38 


.23655996 


10 


1 


33830379 


-0.58575152 


-11 


05109579 


-38 


.88877474 


20 


1 


25987331 


-0.66118370 


-11 


12575192 


-38 


.96332180 


50 


1 


25055050 


-0.67043615 


-11 


13498509 


-38 


.97255257 


100 


1 


24998520 


-0.67100089 


-11 


13554949 


-38 


.97311816 


200 


1 


24991013 


-0.67107586 


-11 


13562419 


-38 


.97319073 


300 


1 


24990207 


-0.67108401 


-11 


13563272 


-38 


.97319829 


400 


1 


24990001 


-0.67108599 


-11 


13563474 


-38 


.97320208 



Table 8. The complete Boltzmann virial 



5.3 Further comments on the poles 

As we saw, if the discontinuity in the potential at the cut-off a is not zero, and also not 
infinite, then our asymptotic expansions, for the location of the poles which appear in 
the expression of the S-matrix, show that the number of these poles is infinite. 

Alternatively, we can use an elegant argument, of Newton[0] and Nussenzveig[|]], 
which argues that the function G(k) = F(k)F(—k), involving the Jost function F(k), is 
an entire function of the variable k 2 , of order 1/2, and therefore has an infinite number 
of zeroes, which then again leads to the conclusion that the S-matrix has an infinite 
number of poles. The basic element in both approaches is an analysis of the behaviour 
of the function G(k) for large values of \k\. 

As we can see directly, from our asymptotic expressions, the zeroes of F(k), for 
these large values of |fc|, depend on the value of the potential at a, rather than on the 
values for r from to a. Thus, as already noted by Nussenzveig, the resulting poles do 
not have much physical significance. He shows that a Yukawa-type potential yields a 
branch cut, which if the potential is cut off, is replaced by an infinite set of poles. 

One further point. We can readily understand that some of the poles represent 
bound states and resonances. In general, it is subtle to understand the physical sig- 
nificance of the poles. Nussenzveig dedicates a chapter in his book (Causality and 
Dispersion Relations) to this purpose. We commend it to our readers. 

6 Conclusion 

We think that the development of a formalism for the second virial coefficient, in terms 
of the poles of the S-matrix, is an attractive one. It is a formalism equally as powerful 
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as the more conventional one, based on phase shifts and bound states, but, for example, 
treats the phase shifts and bound state contributions in a unified way. 

For hard spheres, the number of poles for each i is finite and we see that, with a 
handful of these poles, we can reproduce results obtained by previous methods. 

For other more realistic potentials, and our hard sphere + square well, the number 
of poles is infinite and the phase shifts and the virial converge slowly in terms of the 
poles. We have, however, been able to devise tricks to accelerate this convergence. We 
have also perfected and extended the use of asymptotic expressions for the location of 
the poles. The result is that we can still obtain (and have obtained) results with a 
modest (non forbidding) number of poles. 

These days, with the abundance of numerical power available, and the new methods 
that have been developed to locate poles for realistic potentials, we are freer to choose 
the methods that we might use for virial calculations. 
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Appendix A 

The integrals in (5) are of the type 



f 

Jo 



dk e~ 2(3k 



i 



'it 



£■2 _ u ' 



where we drop, for convenience, the subscript £ in the expression of the k n /s. 
Now,let 

/(A) = e-"~ I dk — 



V.2 



The function /(A) satisfies 



^ T/\\ Xk 2 

-/(A) = - e Aft " 



d\ 2^ ' 

We have therefore 

iw = m - ^ [ 

Introducing 

2 f Z 2 

erf(z) = — ^= / e _t) (if , 
V 7T Jo 

we have 

/(A) = fdk^, - ^-erfHv^) 
By using the residues, note that 



2k n 
ni 
2k n 



and thus 

with the sign - when Q(k nj i) is positive and the sign + when $$(k n! i) is negative. 
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Appendix B 



Here, we derive an asymptotic expression for the location of the poles of the S matrix,or, 
equivalently of the zeros of the Jost function (the zeros of the Jost function F(k) are 
the poles of the S-matrix F(—k)/F(k)). We solve a Volterra equation, which generates 
the Jost solution fi(k,r) and, then, the Jost function. The latter Fe(k) = fe(k,a) 
for a potential which includes a hard core and Fe(k) = \im r _>o(—kryf£(k,r)/(2£— 1)!! 
otherwise. We recall that fi(k,r) is defined by the Cauchy condition: 



t 



lim fi(k,r) exp(-ikr) = i 

This gives us a Jost function which is analytic in the upper half plane > 0. 
For i = 0, let be g(k, r) = fo(k, r) exp(-ikr) . The function g satisfies 

j poo . . 

9(k, r) = 1 - — jf (e**< r ^ - l) V(r') g(k, r') dr' . (fil) 

As usual]?], writing g(k, r) as a series r) = J29n(k, r), we note that it is abso- 
lutely and uniformly convergent when the potential V satisfies / °° r V(r) dr < oo. We 
construct the <? n by the recursive procedure 

g (k,r) = 1 

9n(k,r) = /°° (e 2ife(r '~ r) - l) V(rO Sn-i(*, r') dr' n > 1 . 

We then obtain for a finite range potential, infinitely different iable at the left of its 
cutoff denoted a, the terms 



4P \ 2iA; (2ifc) J 

1 ( V{r) V '(r) , V ff (r) 



4/c 2 V w 2i£; (2iA;) s 



and 



2 



1 m«-r) (v(af-3 V ' {( *\^ a) + 



8k 4 V ( 2ik ) 

1 <V(r) rv { r') d r' + 2V2 ^- V '^'" V ^ dr ' + 



8ik 3 \ y J Jr 2ik 



16fc 4 \ w 2ifc / \ 2 ^ 
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+ 



8fc 3 
1 



V{r') \V{r')- 
V'(a 



16A; 4 



V(a) 



2ik 



+ 



V'(r r 
2ik 
V"{a) 



dr' 



2ik 



etc.. and 



32A; 4 



2ik(a—r) 



V(a) (f°V(r') dr'^j 2 + 0(k~ l ) 



(2ikf 



+ 0(Ar 3 ) . 



To obtain the g„'s, we have used partial integration on the factor containing the expo- 
nential and differentiated the term containing the potential. 

If one approximates the function g(k, r) by the sum s(k, r) = g + g\ + g 2 + #3, the 
Jost function is approximated by s(k, 0), for a potential which is finite at the origin, 
and by s(k, a) exp(ika) for a potential which includes a hard core component. We 
obtain, setting M = J° \/(r')rfr', 



1 + 



+ 



4A; 2 



1 (Via) 



2k 



2M- 



V{a), 



VW -2V {a) + 2M V M-2Mi 



Ak 2 \ V(a) V(a) 2 v ' V(a) 
as the condition which will yield the poles of the S matrix, provided that V(a) 7^ 0. 
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Appendix C 



For higher waves we have to deal with the following free Jost solution 

/7T7" , T (i) 



we(kr) = -kr H^ 1/2 (kr) , (CI) 

where Hffi is the Hankel function of the first kind of order v. We have 

Xmxikr) 1 w e {kr) = (-Y (2£-l)\\ , 

r— >0 

and 

lim exp(— ikr) We(kr) = i e . 
The wis are simply given by || 

wt{kr) = i e e ikr P e (kr) , (C2) 
where the P/s denote the polynomial part of the Hankel, i.e. 

We proceed in a manner similar to that used for I = 0. We introduce 

g t (k,r) = f e (k,r)/w e (kr) , (C4) 

where fe(k, r) is the Jost solution having the appropriate behaviour for r tending to 
infinity 

lim i e exp(-ikr) fe(k,r) = 1 . 
Note that the wVs never vanish for k real. The function gg then satisfies 

- P^At') Pt(-kr'))V(r')gi(k, r') dr' . 

We apply the previous procedure (see Appendix B) which consists in using partial inte- 
gration for the factor containing the exponential and differentiating the term containing 
the potential V(r') multiplied by Pf(kr'). 

When the potential includes an hard core component the factors (kr')~ m m > 0, 
occurring in the polynomial Pg, are bounded by (ka)~ m and therefore goes to zero 
when \k\ tends to infinity. 

We then obtain a formula similar to (15) but where successive derivatives of Pi 
appear. This implies additional ^-dependent terms in the expansion in powers of 1/k. 
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For example, for the potential used before (hard core plus square well) we found, 
Ak 2 (a-a) 2 



D 2ik(a—a) 



A 2 



1 + 



x 



k(a-a) U(6 + i; 



A 



+ 
+ 
+ 



1 



2k 2 (a-a) 2 

i 

6k 3 (a — cr) 3 

x 

26 3 (1 + 6) 3 



-A* + 2A< 



x 



6(6 + 1) 



+ A—x 



2b 2 + x 
b 2 (6+ l) 2 
A 2 



2 6(6+1) 6 2 (6+l) 2 

-6(1 + 6) 3 + 186 3 + x(l + 36 - 9b 2 ) - 2x 2 ) 



{2b 2 + x) 



where A = 2-n[a - cr)/A*; x = £(£ + 1) and 6 = a/ {a - a). 

When the potential has no hard core, the factors l/(kr) m m > occurring in the 
Pi's are no longer bounded and we have to reason differently. 

In fact, when r tend to zero the term 

Pe(-kr) 
Pe(kr) ' 

tends to (— ) and, in so far as the leading term of gj> is concerned, we are left with 



9e{k,r) 



2k 



((-)V*( r '- r 'P f (fcr') 



, 2 1 + ikr + . . . 



1 — ikr + . . . 



- Pe(kr') Pt{-kr'))V(r') dr' . 
When \k\ is large, only the behaviour of the potential at its cutoff a dominates 



9t{k,r) 



: V ( a ) 2ika 
Ak 2 



and the leading asymptotic expression is given by solving 



„2ika 



4k 2 



(C5) 



(C6) 



We then recover the alternating sign, which depends on whether I is even or odd, 
mentioned earlier by Nussenzveigjf]. This dependence disappears when the potential 
incorporates a hard core. 
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Appendix D 



We examine the low temperature expansion of the virial for a pure hard core. We start 
from the equation (10), which for the hard sphere reads: 



[B2) Boltz = -2 l ' 2 \%U £(2£+l) 



7fc + ^? exp( "^^ )erfc( ^ M 



(Dl) 



The equation (Dl) involves the function exp(— z 2 ) eric(iz), where z = Xt k n ^/y/2ir. 
This latter has the asymptotic expression for \ T /a (or equivalently z ) large 



exp(— z 2 ) erfc(iz) 



i + E 



(2j - 1)!! 



(D2) 



Incorporating (D2), written for z = Xt k nt e/y/2Tr, into (Dl) we have: 



Boltz 



-2^ 2 X%M 



E 



V2X T n=1 



V2X T l^t ^ ( V2Xq 

J_ " (2j - 1)!! 

k o ^ h 2j+l \ 2j 

K>n,t j=i K nt A T 



In the previous equation use is made of the property for hard spheres the S matrix has 
no poles for I = and exactly £ poles for I ^ 0. 

The expression for the virial, divided by its classical limit, as in (30) reads: 



(B* 2 ) 



Boltz 



_3_ Xt 
2tt \~o~ 



i + £(2£ + i) i + Ettt^ 

a \ 2 i 



f rr \ l 3 

+ £(2j-l)!!(-)'V - £ 

7=1 KA ? / n=l 



2J+1 



,(^n,£ 0") 



The calculation of the virial requires the knowledge of the sums 



•Sj - ,* = E 



1 V Cr) , 



J<1 



(D3) 



The poles fc n ^ of the S matrix are the zeros of the polynomial part of the Hankel 
function (Eq.(C3) for r = a). Introducing x n ^ = ij{k n ,t cr), these latter are roots of 
the polynomial 



Clm. J 



m=0 



with 



p(x) = E 

(I + m 



dm/ 



2 m (I - m)\m\ 



{DA) 
(£>5) 
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The sums Sj t e, Eq.(D3), are given by 



3<1, 



n=l 



in terms of the roots of the polynomial Eqs. (D4,D5). 
They obey the recursion formula 



Si, 



3,t 



3-1 

-E 



a e- m ,e 



S m ,e — j 



Qi-i,l 
ae,e 

a>e-j,e 



m =i a e,e a e,< 
From (D5) and (D7) we find, after calculation, 



J>2 



I 




Sty 


Sty 


Sj,t 


Sty 


S\ty 


1 


-1 


-1 


-1 


-1 


-1 


-1 


2 


-1 





1/9 


1/27 





-1/243 


3 


-1 








-1/225 


-1/1125 


-1/16875 


4 


-1 











1/11025 


1/77175 


5 


-1 














-1/893025 



The sums S 2 j+ty,j ^ are zero for £> j + 1. 
Rewriting (B£) Boltz in terms of the Sj/s 



(B 



2) Bolt z 



3 A?n 
2^ ^ 



l + £(2j-l)!! 

3=1 



a \ 2 i 



£^0 



and, taking into account the results depicted in the table, we have 



(£>6) 



(D7) 



, (DS) 



Boltz 



+ 



3 

2^ ^ 
472102277 



l + 3vr 



71 



a 
Ay 

(7 



22 9 
3 VA 



1921 
^5~ 



-7T 



(j\ 6 165673 
AW 525" 



-7T 



0" 

A7 1 



10 



+ . . . 



165375 " VA T , 
We thus recover the first terms extracted by BLK. 
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